Anharmonicity due to Electron-Phonon Coupling in Magnetite 



Moritz Hoesch, 1 Przemyslaw Piekarz, 2 Alexey Bosak, 3 Mathieu Le Tacon, 4 
Michael Krisch, 3 Andrzej Kozlowski, 5 Andrzej M. Oles, 4 ' 6 and Krzysztof Parlinski 2 

1 Diamond Light Source, Harwell Campus, Didcot 0X11 ODE, Oxfordshire, England 
2 Institute of Nuclear Physics, Polish Academy of Sciences, Radzikowskiego 152, PL-31342 Krakow, Poland 
3 European Synchrotron Radiation Facility, 6 rue Jules Horowitz, F-38043 Grenoble Cedex, France 
4 Max- Planck- Institut fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 
''Faculty of Physics and Applied Computer Science, AGH- University of Science and Technology, 
Aleja Mickiewicza 30, PL-30059 Krakow, Poland 
'Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, PL-30059 Krakow, Poland 

(Dated: April 19, 2013) 

We present the results of inelastic x-ray scattering for magnetite and analyze the energies and 
widths of the phonon modes with different symmetries in a broad range of temperature 125 < T < 
293 K. The phonon modes with X4 and A5 symmetries broaden in a nonlinear way with decreasing T 
when the Verwey transition is approached. It is found that the maxima of phonon widths occur away 
from high-symmetry points which suggests the incommensurate character of critical fluctuations. 
Strong phonon anharmonicity induced by electron-phonon coupling is discovered by a combination 
of these experimental results with ab initio calculations which take into account local Coulomb 
interactions at Fe ions. It (i) explains observed anomalous phonon broadening, and (ii) demonstrates 
that the Verwey transition is a cooperative phenomenon which involves a wide spectrum of phonons 
coupled to the electron charge fluctuations condensing in the low-symmetry phase. 

PACS numbers: 63.20.dd, 63.20.dk, 71.30.+h, 75.25.Dk 



Discovered in ancient Greece, magnetite played a cru- 
cial role in the history of magnetism. Apart from its 
magnetic properties, which found numerous technological 
applications during the last century, magnetite exhibits 
extraordinary behavior connected with the Verwey tran- 
sition P]. The most prominent feature observed upon 
lowering the temperature is the discontinuous reduction 
of the electric conductivity by two orders of magnitude 
at Ty = 122 K. Extensive experimental and theoretical 
studies of magnetite during last decades revealed very 
complex interrelations between the electronic and struc- 
tural degrees of freedom, which both participate in the 
Verwey phase transition [21 [3]. In spite of this progress, 
the transition is not yet fully understood and remains a 
very intriguing phenomenon in condensed matter physics. 

At the Verwey transition, the crystal structure changes 
from the high-temperature cubic FdZm to the low- 
temperature monoclinic Cc symmetry [JJ . It was first be- 
lieved that octahedral B sites, with mixed-valence Fe 2 5+ 
ions turn to distinct Fe 2+ and Fe 3+ ions below Ty. De- 
tailed diffraction studies of the low-temperature phase 
revealed a much more complex charge-order pattern [5], 
with a wide distribution of valences at the octahedral 
sites [6], and an orbitally ordered pattern distributed 
over three Fe sites (trimerons) [7]. These observations 
are consistent with ab initio calculations [5] [5] and res- 
onant x-ray scattering studies [TUrH5] . that additionally 
suggest the charge-orbital (CO) order appearing a few 
Kelvins above Ty [T3] and possibly existing in the dy- 
namic form even at higher temperatures [15] . Since the 
existence of CO order above the Verwey transition has 
been questioned by other studies [T5] and its role in this 



transition is under debate at present [13], experiments 
performed above Ty might resolve this controversy. 

The key element, which stabilizes the CO order and 
drives the crystal symmetry change is the coupling be- 
tween electronic and vibrational degrees of freedom. Such 
cooperative nature of the Verwey transition is supported 
by numerous experimental observations: oxygen isotope 
effect [17], critical softening of the C44 elastic constant 
[T5] , critical diffuse scattering QjJJ [20] , phonon anomalies 
measured by the Brillouin |21j . Raman [25], and nuclear 
inelastic scattering [23] [23] . The instability of the elec- 
tronic structure is intimately connected with a lattice 
deformation, with certain similarity to the Peierls model 
|25j . In the simplest model, one can explain the doubling 
of the unit cell along the c direction by the condensation 
of the acoustic A 5 mode at q = (0, 0, ^) [25], whereas the 
monoclinic distortion requires at least two primary order 
parameters with A5 and A3 symmetries [27; . 

In this Letter, we report a lattice dynamics study of 
magnetite over a wide range of temperatures above Ty 
and address the question whether or not the electron- 
phonon coupling in the presence of CO fluctuations is 
the principal cause of the Verwey transition. Our study 
reveals increasing anomalous broadening of the lowest 
transverse acoustic (TA) modes along the [001] (A) and 
transverse optic (TO) modes along the [110] (£) direc- 
tions as temperature decreases towards Ty. We show 
by density functional theory (DFT) calculations that the 
anharmonic phonon behavior discovered here: (i) is a di- 
rect consequence of the electron-phonon coupling, and 
(ii) occurs only in presence of strong electronic correla- 
tions which is a novel aspect of the Verwey transition. 



2 



□ Q = (4 4 5) data 
• Q = (5 5 4) data 
X. contribution 




(a) 



T 
10 

energy loss (meV) 



FIG. 1. Inelastic x-ray scattering spectra of FesCU at T = 150 
K. At the X-point (554) scattering from the X4 phonon and 
at (445) from the X3 phonon with a small contribution from 
X4 are dominant. Energy resolution is 1.6 meV. 



Momentum resolved vibrational spectra were obtained 
by inelastic x-ray scattering (IXS) at beamline ID28 at 
the European Synchrotron Radiation Facility (ESRF) 
|28j . Two settings with energy resolution of SE =1.6 
meV and 3.2 meV were used. The sample was a single 
crystal of magnetite, skull melter grown at Purdue Uni- 
versity and annealed for stoichiometry [55]. The sample 
temperature was measured using a Si diode sensor on the 
copper cold finger of the closed cycle cryostat. 

The momentum points Q to be measured were cho- 
sen from the calculations presented below. From the 
calculated phonon Eigenvectors the IXS intensity was 
predicted and favorable X-points Q were selected from 
the highest contrast against other modes and high to- 
tal intensity. The most suitable X4 mode was found 
at Q = (554). The optimal X3 mode was found at 
Q — (445) with a slight contribution from the X4 mode in 
the spectra. The lowest optical A5 mode was measured 
along the line from Q = (333) to Q — (334). 

Raw data with SE — 1.6 meV at the Brillouin zone 
boundary point X are shown in Fig. [l] The lowest en- 
ergy X4 phonon displays a significant broadening, while 
the lowest X3 phonon is only slightly broader than the 
experimental resolution. As the temperature is lowered 
towards Ty the diffuse scattering intensity at zero energy 
loss strongly increases. The phonon energies on the other 
hand show no readily visible temperature dependence. 

To extract precise phonon energies and peak widths 
the spectra were numerically fitted by a superposition 
of peaks derived from the measured resolution function. 
For visibly enlarged phonon peaks these were convoluted 
with a Lorentzian function. For the quasi-elastic inten- 
sity no convolution was applied. The spectra at (554) 
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FIG. 2. (a,b) Phonon energies determined by peak fitting as 
a function of reduced wave vector for branches connected to 
Xa, measured in the vicinity of (554), for two temperatures 
T = 293 K and T = 129 K. Calculated phonon dispersions 
are plotted with solid lines. (c,d) Phonon widths of branches 
connected to X4 as a function of reduced wave vector. 



of Fig. [T] were thus analyzed first, and the determined 
phonon energy and width were kept fixed for the fit of 
the data at (445). The free parameters were the intensi- 
ties and one phonon energy and width per spectrum. We 
found that the broadening had to be applied to the X4 
mode and all connected branches, while the X3 mode and 
connected branches remained close to resolution limited. 

The thus determined phonon dispersions from data at 
SE = 3.2 meV are shown in Figs.[2^a,b) for branches con- 
nected to X4 and in Fig. [3] for X 3 , Fig. [3] also includes 
room temperature inelastic neutron scattering data |30j 
and the acoustic E3-TA mode is fit by a model sinu- 
soidal dispersion, which reveals a downward bending of 
the dispersion by less than 1 meV about half way be- 
tween K and T. The dispersions match within the error 
bars both at room temperature and T = 129 K, with 
only a slight hardening at low temperature visible for 
the X4 mode, in agreement with the measured Fe den- 
sity of states [23l [24] (below Ty the multi-twinned nature 
of the sample prevents the determination of the phonon 
dispersion). In the following the A 5 -TA, Ei-TO and E 2 - 
TO will be denoted as the branches connected to X4, 
and A5-TO and E3-TA as connected to X3. In contrast, 
the phonon widths in Figs. [2jc) and [2^d) have a much 
stronger temperature dependence. The width increases 
on lowering the temperature which contradicts the usual 
situation, where anharmonicity increases with increasing 
temperature due to phonon-phonon interactions. 

The full temperature dependence of selected modes, 
measured with SE = 1.6 meV is shown in Fig. [4] Re- 
markably, the modes connected to X4 show an increas- 
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FIG. 3. Phonon energies as a function of reduced wave vector 
for branches connected to X3, measured in the vicinity of 
(445) and (334). Data at two temperatures T = 293 K and 
T = 129 K are shown. Calculated phonon dispersions are 
plotted with solid lines. Data from neutron scattering at room 
temperature [3U] and a sinusodal fit to represent a typical 
phonon dispersion are also shown. 



ing width on lowering the temperature towards Ty. The 
X4 and A 5 -TA phonons show already anomalous widths 
at room temperature moderately increasing with decreas- 
ing temperature. After reaching maximum values around 
150 K, their widths are reduced again in a temperature 
range of about 20 K above Ty. The Si and £2 modes, 
measured at the K point, show even stronger and linear 
increase of widths over the entire range of temperatures. 

The DFT calculations were performed using the VASP 
program |31j within the generalized gradient approxima- 
tion (GGA) and the full-potential projector-augmented 
wave method [32]. Phonon energies for the cubic Fd3m 
symmetry were calculated in the lxlxl supercell with 
56 atoms using the direct method [33], implemented in 
the PHONON program [34] . On the basis of the Hcllmann- 
Feynmann theorem, the atomic forces were obtained by 
displacing atoms from their equilibrium positions. For 
the cubic symmetry only three independent displace- 
ments of Fe(A), Fe(-B), and O atoms are sufficient to 
derive all force-constants matrix elements and the dy- 
namical matrix. In order to study anharmonic effects, 
we have calculated the Hellmann-Feynmann forces and 
phonon energies for two independent sets of displace- 
ments with the amplitudes u = 0.02 and 0.04 A. The 
effect of local electron interactions on phonon energies 
was investigated by performing calculations within the 
GGA+C7 approach, where the on-site Coulomb interac- 
tion and Hund's exchange parameters, U = 4.0 eV and 
J = 0.8 eV, are the same as in the previous studies [2"7] . 

Considering the A and E directions, we found a very 
good agreement for both the A 5 modes between the ex- 
perimental points and the calculated phonon dispersions 
for u — 0.02 A and only a small shift of the Si mode 
to lower energies (Fig. [2]). The largest discrepancy is 
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FIG. 4. (a,c) Energies of phonon modes at the -ft'-point as 
a function of temperature: (a) A5, X4 and X3, and (c) £1 
and £2; Widths of these phonon modes are displayed in (b) 
and (d), respectively. Straight dashed lines in (b) and (d) 
were fitted to points at high temperature. Widths within the 
hatched area are not reliably extracted from the data. 



found for the S2 mode, which shows anomalous softening 
induced by the local electron interactions [27], directly 
connected with the anharmonic behavior discussed be- 
low. The S3 mode exibits anomalous downward bending, 
much larger than those observed in the experiment. 

Coulomb interactions strongly influence phonon dis- 
persions. For {7 = [Fig. [5]fa)] the phonon energies ob- 
tained for u = 0.02 and 0.04 A are very close to each other 
and differ by less than 0.5 meV (not shown), which means 
slight hardening of phonons for increasing displacement 
amplitudes. It demonstrates that the system is harmonic 
and phonon frequencies depend only very weakly on vi- 
brational amplitudes. The situation changes in a dra- 
matic way in presence of local electron interactions [see 
Figs. p)Jb,c)] — for U = 4 eV the phonon energies depend 
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FIG. 5. Calculated phonon dispersions of FesCU obtained 
from GGA+U combined with the phonon software [33], for: 
(a) (7 = 0, and (b) (7 = 4 eV, shown by solid (dashed) lines 
for a displacement of it = 0.02 A (0.04 A). Selected branches 
connected to the lowest X3 and X4 optical modes are high- 
lighted. The difference between the phonon dispersions ob- 
tained at both values of u for (7 = 4 eV is shown in (c). 



strongly on u and substantially soften for u = 0.04 A. The 
largest changes are observed for the branches connected 
to AT 4 with a broad maximum in the energy differences 
at the X point. A weaker but still pronounced effect is 
visible also for branches connected to X3. 

We emphasize that the same modes which have anoma- 
lous temperature dependence of phonon widths exhibit 
strong anharmonic behavior in the theory when (7 = 4 
eV. Moreover, the phonon widths for the A5-TA and £2- 
TO modes (Fig. [2]) show similar trends as the plot of en- 
ergy differences, see Fig. [5^c). In both cases the largest 
values are located in the vicinity of the X point. Along 
the A direction, one finds a pronounced increase of val- 
ues in both figures around q — (0,0, \). The phonon 
widths at room temperature have a maximum around 
q = (0, 0, 0.7) and decrease systematically when q in- 
creases towards the X point, and further along the £ di- 



rection. At low temperature a two-peak structure for the 
(^-dependent width is apparent: one peak at q ~ (0, 0, 0.8) 
and the other one at the K point. 

The observed temperature dependence of the phonon 
widths and the exact positions of their maxima cannot 
be properly described within the present calculations, see 
Fig. [5jc) . The shift of peaks away from the X point 
indicates the existence of incommensurate correlations, 
which are missing in the supercell calculations due to the 
geometrical constraints. 

In standard theories, large phonon widths indicate ei- 
ther an anharmonic behavior or strong electron-phonon 
coupling [35] . Our ab initio study demonstrates that 
these two effects are intimately connected in magnetite 
and both result from local electron interactions. These 
interactions induce polarization of the minority-spin 
partly occupied t2 g orbitals. In such a correlated state, 
the coupling between electrons and lattice enhances and 
stabilizes the CO order, effectively reducing the total en- 
ergy \27\ I36j . This in turn modifies the interatomic inter- 
actions and generates the anharmonic potential. Indeed, 
an anharmonic double-well potential has been obtained 
for the A5 and X3 modes in the presence of strong elec- 
tron correlations [57]. This mechanism leads to the Ver- 
wey transition with the static structural distortion and 
frozen CO order below Ty. In contrast, above Ty the 
CO degrees of freedom are dynamic and couple strongly 
to phonons reducing their lifetimes. Therefore we con- 
clude that the anharmonicity in magnetite is induced by 
the coupling between phonons and CO fluctuations which 
leads to phonon broadening [37] . 

We remark that a strong electron-phonon coupling is 
consistent with the polaronic short-range order observed 
above Ty [38-440] . The anomalous phonon broadening 
correlates also very well with the critical diffuse scatter- 
ing existing over a broad region in reciprocal space and 
over a broad temperature range [20 . In both cases, the 
strongest response is observed at incommensurate wave- 
lengths in the reciprocal space. Charge fluctuations ob- 
served in a wider temperature range may be responsible 
for distinct properties of the Verwey transition in non- 
stoichiometric or doped magnetite (24] SU [42] . 

In summary, we have presented the IXS measurements 
of the lowest phonon modes along the A and £ direc- 
tions in magnetite above the Verwey transition — they 
reveal an anomalous behavior of phonon widths that in- 
crease with lowering temperature. Their maxima are 
found away from the high-symmetry points due to the in- 
commensurate character of critical fluctuations. Within 
the ab initio study we have found that the same phonon 
modes show strong softening for larger amplitudes in 
the presence of local electron interactions in the Fe(3d) 
states. This phonon anharmonicity is induced by the 
charge-orbital fluctuations above Ty — it demonstrates 
that the Verwey transition is triggered by a cooperative 
electron-phonon mechanism at large U [27j . 
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